Entangling two distinguishable matter-wave bright solitons via collisions 
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We investigate numerically the collisions of two distinguishable quantum matter- wave bright soli- 
tons in a one-dimensional harmonic trap. We show that such collisions can be used to generate 
mesoscopic Bell states which can reliably be distinguished from statistical mixtures. Calculation 
of the relevant s-wave scattering lengths reveals that such states could potentially be realized in 
quantum-degenerate mixtures of 85 Rb and 133 Cs. In addition to fully quantum simulations for two 
distinguishable two-particle solitons, we use a mean-field description supplemented by a stochastic 
treatment of quantum fluctuations in the soliton's center of mass: We demonstrate the validity 
of this approach by comparison to an effective potential treatment of the quantum many-particle 
problem. 
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Generating quantum entanglement between meso- 
scopic objects over mesoscopic distances allows explo- 
ration of a fascinating "middle-ground" between quan- 
tum and classical physics [U [2] ; such mesoscopic quan- 
tum superpositions are relevant for quantum-enhanced 
interferometry [3J. Quantum-degenerate gases offer a 
highly controllable platform for such studies: Scattering 
quantum bright matter-wave solitons fl"lU2] in quasi-one- 
dimensional (ID) trapping geometries has been suggested 
as a way to create mesoscopic entangled states in single 
species Bose-Einstein condensates (BECs) [T3J Q3] . 

In this Letter we consider a dual-species BEC Q~51 [TB] , 
and show that collisions of distinguishable quantum 
bright matter-wave solitons can be used to generate 
mesoscopic Bell states [17] (cf. [18]). 
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where \A,B) (\B,A)) signifies that the BEC A is on the 
left (right) and the BEC B is on the right (left). In par- 
ticular, we show that a favorable combination of inter- 
and intra-species s-wave scattering lengths means that 
such states may be realized using 85 Rb and 133 Cs mix- 
tures. We also show that the interference properties of 
these bright-soliton Bell states distinguish them from sta- 
tistical mixtures. In contrast to the Bell ground states 
associated with double-well potentials, our collisionally- 
generated Bell states are robust to the presence of asym- 
metries. While distinguishable solitons are essential to 
produce Bell states, entanglement generation for solitons 
of the same species was investigated in [3]. 



For our quasi-lD system, we consider an experimen- 
tally motivated harmonic confinement lu = 2?tf. Mix- 
tures of ultracold gases can be confined in a common 
optical trap with the same trap frequencies |19j , yielding 
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where ttia (tub) is the mass of an atom of species A (B); 
the interactions g are related to the scattering lengths a 
and the perpendicular trapping- frequency, fx, via [20] 



= hfxa . 
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We use the Lieb-Liniger model [21] for two species with 
additional harmonic confinement 
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where <?a < (c?b < 0) is the intra-species interaction of 
species A (B), ^ab > the inter-species interaction and 
Xj (yj) are the coordinates of atoms of species A (B). 
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We suggest to prepare the two solitons independently; 
for weak harmonic confinement a single soliton has the 
ground state energy (cf. [22] ) 



1); Se{A,B}. (5) 
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Thus, our system has the total ground-state energy 

E = E A (N A ) + E B (N B ) . (6) 

The total kinetic energy related to the center-of-mass 
momenta HK$ ( S G {A, B}) of the two solitons reads 



kin 2N A m A + 2N B m B ' ^ 

We extend the low-energy regime investigated for 
single-species solitons in Refs. [7J [TSl [53] to two species: 

E kin < min{A A , A B }, 
A s = \E S (N S - 1) - E S (N S )\ . 

In this energy regime, each of the quantum matter-wave 
bright solitons is energetically forbidden to break up into 
two or more parts. Highly entangled states are charac- 
terized by a roughly 50:50 chance of finding the soliton A 
(B) on the left /right combined with a left/right correla- 
tion close to one indicating that whenever soliton A is on 
the one side, soliton B is on the other: 



T («5) 



dxi 



d.r 



N A 



dyi 



d.r 



n a 



OG 
OO 



dyi 



dy NB \*\ 2 

dy NB \y\ 2 , 
(8) 



where = $(xi, . . . , xjv a , y\, . . . , Vnb) ^ s the many- 
particle wave function (normalized to one) and S > 0. 
The correlation ^(5) will serve as an indication of en- 
tanglement: Bell states |l]) are characterized by 7 ~ 1 
combined with a 50:50 chance to find soliton A either on 
one side or on the other. 

We begin by investigating entanglement-generating 
collisions of two distinguishable two-particle solitons 
(dimcrs). Discarding cases where the two solitons have 
distinct total masses N A m A and N B m B (small differences 
in the total masses would introduce small asymmetries 
without changing the physics), leads to m A = m B = to, 
which corresponds to having two hyperfine states of the 
same species. To describe the collisions of the two dimers, 
we discretize the Hamiltonian Q, yielding the Bose- 
Hubbard Hamiltonian (cf. [2"4"] 1 
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FIG. 1. (Color online) Collisions of two distinguishable dimers 
in the Bose-Hubbard Hamiltonian (TsJ) . (a) Single particle 
density g(x) of dimer A in a two-dimensional projection as a 
function of space and time (r is the oscillation period without 
inter-species interaction, Ua ~ — 3J, Ub = — 3J, Uab ~ J 
and C = 0.002 J), (b) Single particle density of dimer B, 
parameters as in panel (a), (c) The same dimer as in panel (a) 
but the wave function is numerically turned into a statistical 
mixture at t — r/2. (d) Center-of-mass density qcom(X — Y) 
if the inter-species interaction is switched off at t — r/2, all 
other parameters as in (a) and (b). 



where U Al U B and U AB are the intra-species and inter- 
species interactions, the hopping is given by J ~ 
H 2 /(2md 2 ) for grid spacing d — > and C = 0.5tou; 2 gP . 

Figure [T] shows two-dimensional projections of the dy- 
namics of two distinguishable dimers. The two dimers 
were numerically prepared in the ground state of two spa- 
tially separated harmonic oscillators via imaginary time- 
evolution [25]. At time t — they were transferred into 
the same harmonic oscillator potential (without over- 
lap). Subsequently, the time-evolution was calculated 
using the full Schrodinger equation corresponding to the 
Hamiltonian After the first collision, a measurement 
would reveal dimer A on the left and dimer B on the right 
or vice versa [the correlation ^ is j(d/2) ~ 0.988]. 

As the sizes of the dimers in panels Fig. [T] (a) and (b) 
are not too large compared to the oscillator length, after 
the second collision both dimers are more likely to be on 
the side opposite to their initial condition than at the 
same side (cf. the single soliton case 23J). This can be 
used to distinguish a pure quantum superposition from a 
statistical mixture [Fig. [I] (c)]. A more general approach 
extends the center-of-mass density [26] to two solitons: 
first one measures the center of mass X and Y of solitons 
A and B and then plots the resulting density grjoM as 
a function of the difference X — Y . This works both 
for superpositions of plane waves exp[ii^A"] exp[— iKY] + 
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FIG. 2. (Color online) The s-wave scattering lengths for the 
ground state of 85 Rb 133 Cs, 85 Rb and 133 Cs respectively, (a) 
Scattering lengths are calculated using a coupled-channels 
method |16j with a fully decoupled basis set at a collision 
energy of 1 pK. The calculations are performed using the 
MOLSCAT program 27 adapted to handle collisions in ex- 
ternal fields 128] . The RbCs potential is from |16j . Rb from 
[29] and Cs from [30]. Resonances for 85 Rb 133 Cs are at 3.10G, 
4.27G and 6.76 G (b) Zoom of (a). 



exp[— iKX] exp[iK Y] with 

q C om(X-Y)<x{cos{K(X-Y)}} 2 (10) 

and when the two wave packets recombine [Fig. [I] (d)] . 

To show that attractive intra-species interactions and 
repulsive, tunable inter-species interactions are exper- 
imentally feasible, we calculate the s-wave scattering 
lengths for 85 Rb 133 Cs. The results displayed in Fig [2] 
shows a candidate inter-species Feshbach resonance at 
6.76 G suitable for our requirements. For lower mag- 
netic fields the magnetic field can be stabilized to up to 
100 [iG [3T]; shielding allows stabilization to lmG below 
10 G. Although the masses of the atoms A and B now 
differ, we can still have two solitons of roughly the same 
total masses N A m A and N-qtu-q as in Fig. [T] 

Behavior for larger particle numbers can be described 
by the Gross-Pitacvskii equation (GPE) (cf. 32 ) 
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where the single particle density | cps (a;, | 2 is normal- 
ized to Ng (S € {A,B}). Before we reconsider the state- 
ment that the nonlinearity of the GPE does not allow 
to describe quantum superpositions, we investigate the 
GPE-dynamics in Fig [3} For high kinetic energies, both 



solitons split such that roughly 50% of the particles of 
each species are on each side of the center of the trap for 
times near t = T/2; as expected for a GPE-solution, the 
correlation function depicted in panel Fig [3] (c) shows 
no signs of Bell states. For low kinetic energies, sin- 
gle species solitons show jumps in the reflection coeffi- 
cient [T21[231|33] which we have conjectured to correspond 
to quantum superpositions on the A^-particle quantum 
level [23]. As a starting point for investigations of Bell 
states, we thus also look for parameters for which there 
are jumps in the reflection coefficient. We extend the 
exact lower bounds on the reflection coefficient R for a 
single soliton scattered off a barrier [23 to the case of 
two distinguishable solitons [Fig. [3] (d) and (e)], where 



min{i? > 0.5} . 



(11) 



Panels (d) and (e) of Fig. [3] show why, as in Fig. [I] we 
should use two solitons of roughly equal total masses 
N A m A and Nbttib, A^b ~ l-6Ac s for entanglement gen- 
eration: only at approximately this ratio jumps occur in 
both reflection coefficients at reasonable kinetic energies 
corresponding to not too long timescales (without a trap) 
or not too small initial displacements (in the presence of 
a trap). 

To test the conjecture that for jumps in the reflection 
coefficient on the GPE level indeed correspond to quan- 
tum superpositions on the A^-particle quantum level [23j , 
we use the truncated- Wigner approximation, which de- 
scribes quantum systems by averaging over realizations 
of an appropriate classical field equation (in this case, the 
GPE) with initial noise appropriate to either finite [34] 
or zero temperatures [11]. While the GPE assumes both 
position and momentum are well defined, this is not true 
for a quantum particle of finite mass. Our truncated- 
Wigner calculations for the center-of-mass wave function 
use Gaussian probability distributions for both (satisfy- 
ing minimal uncertainty), 

In order to demonstrate that the truncated- Wigner ap- 
proximation is indeed a valid approach to describe the 
short-time behavior of mesoscopic quantum superposi- 
tions, Fig.|4]starts with the case where a light soliton hits 
a heavy, non-moving soliton. In panel (a), the rigorously 
proved [T3J effective potential approach [7] 113^] demon- 
strates the emergence of a Schrodinger-cat state when 



1 This approach replaces the Af-particle Schrodinger equation 
by a single-particle Schrodinger equation for the center-of- 
mass coordinate. The effective potential is the convolu- 
tion of the soliton with the potential seen by single parti- 
cles [7] 1131 . in our case the soliton B which is chosen to 
be a factor of N^g^/(N\g\) = 10 narrower than soliton 
A: [/i[cosh(AT B gBz)]" 2 - 2C/i/(7V B g B )<5(^)- If soliton A 
hits the narrow soliton B [Fig [4] (a)], the effective potential 
C/o[cosh(A r A9A2;)] _2 has the GPE-shape of the A-soliton with 
U = N A UiNj±g A / (A^b^b) = 10C/i, for N A = 100. 
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FIG. 3. (Color online) (a) Single particle density for a GPE- 
soliton of 133 Cs atoms splits at high kinetic energies when 
hitting the soliton made of 85 Rb atoms [panel (b)]. Parame- 
ters are oc s — 982. 5ao, <iRb = — 459.2ao, / = 1Hz, f± — 16Hz, 
Nq s — 17, A^ab — 17, aRbCs = 8351ao, and the starting dis- 
placement from the trap center: ±55.2/im (c) As expected for 
the nonlinear GPE, at t ~ T/2 we have a low corr elation q8p . 
(d) and (e) show lower bounds on 7? m j n [Eq. (Ill cf. [23]" 



the GPE predicts the stepwise behavior of the reflection 
coefficient explained in Refs. QH 123 [33J [Fig|I](c)]. 

In panel (b), the truncated- Wigner approximation 
transfers the analytic approximation for the classical- 
particle-like behavior of the GPE-soliton [35] into a good 
qualitative agreement with the A^-particle predictions in 
panel (a) up to the time where both parts of the wave 
function recombine and quantum interference becomes 
important. 

In order to find parameters for which Bell-states can 
occur, we have to go to lower kinetic energies than in 
Fig. [3] (a) and (b) by choosing a smaller initial displace- 
ment; as suggested by Fig [3] (d) and (e) we further- 
more use A/RhTORh ~ NcstriCs- Applying the truncated- 
Wigner approximation for the center-of-mass wave func- 
tions of both solitons leads to the single-particle densities 
displayed in Fig. |4](e) and (f); due to the low kinetic en- 
ergies the single-particle densities near t « 0.6T should 
be a sign of entanglement. The value of the correlation 
function close to one [Fig. [4] (d)] shows that we indeed 
have found a Bell state. 



FIG. 4. (Color online) (a) Single particle density for an N- 
particle quantum bright soliton hitting a narrow, heavy non- 
moving soliton, computed using the effective potential ap- 
proach (footnote [l]). (b) As in (a) but using the truncated- 
Wigner approximation for the center of mass. The parameters 
are in the low kinetic energy regime such that the GPE pre- 
dicts a stepwise behavior of reflection coefficient as a function 
of the initial displacement shown in panel (c); Uo — 12hio (cf. 
footnote[T]). (d)-(f) Applying the truncated- Wigner approxi- 
mation for the center of mass in the low-kinetic-energy regime 
to the two-species GPE, leads to a much higher correlation Q 
than in Fig. [3](c), thus, indicating a Bell state. Parameters 
as in Fig. |3j except for f± = 70Hz, N Cb ~ 12, N Rh ~ 19 
(corresponding to ATrtbTOrtb ~ Nc s mc s ), aRbCs = 63.6a , and 
the initial displacement from trap center: — 8.7fim for 133 Cs 
[panel (e)] and +8.7^m for 85 Rb [panel (f)]. 



To conclude, based on predictions made on the level 
of many-particle quantum calculations (using the Lieb- 
Liniger model) , we demonstrated numerically that meso- 
scopic Bell states can be generated via colliding two 
distinguishable quantum matter- wave bright solitons. 
Matter-wave bright solitons in 85 Rb- 133 Cs mixtures are a 
promising candidate system for experimental realization 
of such Bell states. We suggest to switch off the inter- 
species interaction once the Bell state has formed and to 
then use interference fringes in the combined center-of- 
mass density ( 10 1 to demonstrate that the result is indeed 
a quantum superposition (as shown in Fig. [I]). 
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